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OO A flexible sheet in uniform parallel flow is studied in order to quantify 

its fluid dynamic drag and fluid-elastic stability characteristics. An experi- 
mental campaign is undertaken that involves a cantilevered sheet in air flow 
characterised by Reynolds numbers of order R = 10^ — 10^. The properties 
of the sheet include: constant mass per unit area; small but finite flexural 
^ rigidity; varying aspect ratios from within the range 0.43 < l/L < 1, where 

r]^ L and / denote the length and width respectively; and tension applied at 

^ the trailing edge. The unique aspect of the present work is an investigation 

_ ^ into the influence of in- plane tension on both the fluid drag and fluid-elastic 

stability of the sheet. In the absence of tension, the configuration resem- 
(-H bles a flag and the drag coefficient is observed to decrease with increasing 

Cin aspect ratio and Reynolds number. In the presence of tension, the fluid 

drag is significantly reduced in the region below the critical flow velocity at 
CN which convected wave instabilities appear. This critical flow velocity can 

^ be increased through the moderate application of in-plane tension. Under 

QQ lateral tension, the drag of the sheet is given to good approximation by the 

OO turbulent boundary layer drag law for a flat plate. Once stability is lost, 

however, the drag coefficient increases rapidly with Reynolds number due to 
convected waves travelling over the sheet's surface. 

Keywords: fluid-elastic stability; fluid drag; uniform flow; flexible sheet; 
^P. in-plane tension 
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S 1 Introduction 



Through casual observation of a flag, otherwise known as a flexible sheet, suspended 
in a breeze, the fluid-elastic instability known as flutter is clearly visible. This 



problem flrst gained attention through Rayleigh (1879) who considered a massless 



tension-free surface of inflnite length in an effort to understand flag flutter. During 
a flutter episode, the flag extracts energy from the fluid, resulting in dispersive 
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waves that originate at the leading edge and propagate towards the traihng edge 
with increasing amphtude. This oscillatory motion consequently affects the ffuid 
dynamic drag which is expected to increase once stability is lost and the frontal 
area exposed to the flow increases. Although a flag in air flow is a familiar example 
of a fluid-elastic system, there are many other examples and applications — the 



propulsion of fish (Daniel, 1984), the dynamics of towed underwater cables (Bowl- 



ing, 1988) and sheets, the biological mechanisms behind snoring (Huang, 1995), 



filaments (Zhang et al. 2000 Schouveiler et al. 2005), paper fiutter in processing 
and printing (Watanabe et al. , 2002), drag stabilisation (Auman and Wilks, 2005), 



and general slender structures in axial flow ( Paidoussis , 2004 ) to name but a few — 
that motivate a study on fluid-elastic instabilities of flexible sheets. However, what 
has received limited treatment throughout available literature, is the influence of 
in-plane tension on a flexible sheet in flow. Additional rigidity through tension is 
expected to stabilise the sheet's fluid-elastic response and consequently reduce its 
fluid dynamic drag. These effects are investigated here. 

A classical flutter-type instability of a flexible sheet of small but flnite thickness 
occurs when some threshold of the incident fluid velocity or perturbation of the 
system is breached. The point at which this occurs is termed the critical velocity Uc- 



In contrast to a forced excitation, such as the reattachment of shear layers (Allen 



and Smits, 2001) or though mechanical means (Taneda and Tomonari, 1974), we 



are concerned with a self-excited response. This is often characterised by convected 



waves or so-called limit-cycle oscillations (see Zhu and Peskin 2002 Tang et al. 



2003 Connell and Yue, 2007, among others). The nature of the problem dictates 



that it is fully coupled whereby the properties of the sheet and the surrounding 
fluid are equally important and intimately linked to the sheet's response. Both 
stability and the subsequent magnitude and frequency of the instability depend 
on a competition between the fluid's momentum and the sheet's inertia, and the 
restorative influences of the sheet's flexural rigidity and any in-plane tension — 
which can be supplied by either the boundary layer or by external means. End- 
point boundary conditions also play an important role; from free-free conditions 



that resemble flsh like swimming (Farnell et al. , 2005), to clamped- free restraint 
indicative of a flag suspended in flow (Fitt and Pope, 2001). 

To aid in our physical understanding, the character of the flow surrounding the 
sheet can be divided into two separate regimes based on the aspect ratio ^ = 
l/L, where / and L denote the sheet width and length respectively. For =2/ <^ 1, 
we essentially have a long slender sheet and the flow about it can be considered 



approximately two-dimensional in the cross-flow direction (Lighthill, 1960). Here, 
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flow separation at the trailing edge is of little consequence compared to boundary 
layer induced tension and added mass effects. In this regime, Uc decreases as 
the aspect ratio increases, and some success has been achieved at predicting Uc 



Lemaitre et al. 


2005 


Souilliez et al. 


2006 



for instance). 

On the other hand, when =2/ = 0{1), the characteristics of the flow are far more 
complicated due to the three-dimensionality of the sheet and the importance of flow 
separation and shear layers in the form of an unsteady wake and downwash over 
the cross-flow direction. In a strictly two-dimensional setting, the effect of an un- 
steady wake on a sheet of infinite width has been considered both theoretically (see 
Huang[ 1995 Zhu and Peskin[ 2002 Argentina and Mahadevan 2005 Tang and 



Paidoussis, 2008, among others) and experimentally in the soap film experiments 



of Zhang et al. (2000). The study by Tang and Paidoussis (2008) demonstrates 
that the wake is of little consequence to stability when the sheet is considered 
lon^ However, for a sheet of finite width, it is generally understood that it will 



be theoretically more stable than its two-dimensional counterpart (Eloy et al. 



2007). The coupling of unsteadiness with three-dimensionality has only recently 



been approached numerically (Tang et al. , 2003) with some success at determining 
the critical flow velocity and the sheet's qualitative behaviour after stability is lost. 

The distinctive feature of a flexible sheet is its low flexural rigidity. This feature 
complicates any experimental study where it is important for the sheet to maintain 
a planar form with little or no geometric irregularities as the flow velocity is in- 
creased up to, and beyond, the critical velocity. To circumvent this, most authors 



choose to study a sheet mounted vertically (Taneda, 1968 Datta and Gottenberg 



1975 Lemaitre et al. , 2005) allowing gravity to impose a linearly varying axial 



tension that reduces to zero at the trailing edge. This gravity-induced in-plane 
tension dominates over the minimal amount provided by fluid friction under real- 



istic flow conditions (Lemaitre et al. , 2005). Alternatively, one may mount a sheet 
horizontally in flow and employ sufficient in-plane tension or so-called pre-tension 



(Coene, 1992 Watanabe et al. , 2002 Yamaguchi et al. , 2003 Morris- Thomas and 



Steen , 2008 ) to ensure a planar sheet form up to the critical flow velocity. 



Although applying in-plane tension to a flexible sheet is convenient for experi- 
ment, it does however replicate some important practical applications; for instance, 
print media and paper manufacture, and aircraft towed advertising banners. In 
printing and paper manufacture, knowledge of the stability characteristics of the 



^The term 'long' is subjective in a two-dimensional setting. Tang and Paidoussis (2008 1 cite a 



mass ratio of pL/M > 4 to distinguish a long two-dimensional sheet 
and M the sheet's mass per unit length. 



p being the fluid density 
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sheet are important; for towed banners on the other hand, fluid dynamic drag re- 
duction assumes priority over convected ripple like mode^ In the context of pure 



transverse tension (cross-flow normal) applied to a simply supported sheet, Watan- 



abe et al. (2002) demonstrates an elevated critical velocity by suppressing lateral 



sheet displacements. Moreover, Watanabe and colleagues present some revealing 
illustrations that show the flow surrounding a sheet in the presence of pure lateral 
restraint to be largely two-dimensional, distinctly unseparated, with progressive 
flexural waves exhibiting a remarkable similarity to those of a cantilevered pipe 
conveying fluid (Watanabe et al. , 2002). In contrast, but with similar conclusions, 
Coene (1992) examined a free-free paper sheet of ^ = (9(0.1) with tension sup- 



plied axially (normal to the trailing edge). A clear increase in the critical velocity 
is shown. However, given that the sheet's flexure was not defined, coupled with 
the inherent anisotropy of paper, any direct comparison with existing or new data 
is speculative at best. 

The fluid dynamic drag of a sheet is expected to increase once stability is lost 



to flutter (Taneda, 1968). This is because: the frontal area exposed to the flow 
increases in the presence of flexural waves; and the wake transforms from a well de- 



fined von Karman vortex street to one comprising large eddy formations (Taneda 



1968, especially Fig. 13). While this phenomenon has been demonstrated for lami- 



nar flows with in-plane tension provided by gravity (Taneda, 1968), its observation 
in turbulent flows of more practical interest has suffered because of geometric irreg- 



ularities of the sheet up to the point of instability (see Fairthorne 1930 ; Auman 



and Dahlke 2001; Carruthers and Filippone 2005). Like a flat plate, the drag 



coefficient of a flexible sheet decreases with increasing aspect ratio (Taneda, 1968 



Auman and Dahlke 2001 Carruthers and Filippone, 2005); and, furthermore, we 



can expect its magnitude to approach that of a flat plate under sufficiently small 
elastic deformations. To the authors' knowledge, the effect of additional rigidity 
through in-plane tension, although likely to postpone an expected rise in the drag 
coefficient at the onset of a fluid-elastic instability, has hitherto not been studied. 

The present paper examines a rectangular flexible sheet of small flexural rigidity 
in uniform parallel flow. The leading and trailing edges of the sheet comprise 
clamped and free boundary conditions respectively. Our primary goal is to address 
the issue of in-plane tension and its impact on both the fluid-elastic stability and 
fluid dynamic drag characteristics of the sheet. An experimental campaign is 
undertaken whereby the sheet's dimensions are chosen such that 0.43 < ^ < 1. 



^These ripple like modes resemble the wrinkles which appear across the surface of membranes 
under in-plane tension (see Wong and Pellegrino 2006 1 . 
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In-plane tension is applied to the sheet at a given angle 6 to the trailing edge. 
The magnitude of this tension is varied from zero (a flag) to values in the region 
where tension induced rigidity is expected to dominate over flexure. The instability 
we are concerned with is of the flutter-type which is characterised by convected 
waves propagating from the leading to the trailing edge of the sheet's surface. 
However, under certain conditions transverse mode shapes are possible. The paper 
is organised as follows: appropriate scalings and an approximate analytical model 
describing the instability of a flexible sheet are presented in § |2j the experimental 
campaign undertaken is described in § |3] results pertaining to the flutter-type 
instability § 4.1| the fluid dynamic drag § 4^ and the fluid dynamic drag in the 



presence of in-plane tension § |4.3 are then discussed; and finally, conclusions are 
provided in § |5| 

2 Scaling and a qualitative model 

We first consider a simple two-dimensional analytical model of a flexible sheet. This 
enables us to construct a consistent set of dimensionless parameters, and facilitates 
an understanding of how these parameters affect a flutter-type instability. This 



approach is motivated by Coene (1992) and Shelley et al. (2005) and we combine 
their ideas to include linear elastic deformation of the sheet due to both flexure 
and axial tension. 

We consider a flexible sheet of length L comprising an isotropic material of mass 
per unit area M immersed in an ideal fluid in the presence of uniform parallel flow 
of velocity U . The flow approaches from X = — oo and we assume that the sheet 
occupies the plane Y = in its quiescent position. The vertical displacement of 
the sheet from its mean position is denoted by w{X, t), which we assume small — 
w <^ A, A being a typical wavelength, or w <^ L, L being the sheet length. This 
assumption is justified given that we are only interested in the stability criteria, 
rather than post-stability behaviour, of the sheet. Under the restriction of linear 
elasticity, the vertical displacement of the flexible sheet is governed by 

which is essentially the classical Euler-Bernoulli beam equation with an additional 
term accounting for axial tension per unit width of the sheet T. The fluid pressure 
acting across the sheet is denoted AP{X]t) and B = Eh^ /12{1 — z/^) defines the 
flexural rigidity, where h is the sheet thickness, E the elastic modulus, and u 
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denotes Poisson's ratio. 

We now assume a harmonic response of the sheet, and that the amphtude of 
this response is much smaller than the sheet length. Its linear displacement from 
the mean position Y = can then be written 

w{X,t) = Aexp{iKX -int), (2) 

where K denotes the spatial wavenumber and Q the circular wave frequency. In 
adopting (|2|, we are essentially implying that \ <^ L, which is an obvious sim- 
plification. However, this is justified given our motivation involves identifying 
consistent scaling parameters and deducing their qualitative contribution to the 
sheet's stability. 

Under the assumption that the sheet motion and zero flux boundary condition 
are linear about the plane Y = 0, the appropriate form of the fluid pressure, in 
accordance with ([2]), acting across the sheet is (cf. Lamb, 1932, Art. 232) 



AP(X; t) = 2-^{n - UKf w, (3) 
\K\ 

where = y/lC^. In scaling the governing equation, it is convenient to adopt: 

X = xL, w = 7]L, t = T^ ^=T' fi = ^Y, AP = pU^Ap. (4) 

U Ij Ij 

Introducing these into the governing equation ([T]) with ^ provides 

Vtt + \vxxxx - i^Vxx = '2Tri{i^ - k)'^'n, (5) 

\k\ 

which describes the displacement of the sheet in terms of three dimensionless quan- 
tities defined as follows: 

The first parameter, /x, describes the ratio of fluid forcing to the inertia of the 
sheet with time scaled as t{U/L), and the second, v, is commonly referred to as 
the reduced velocity and is essentially the ratio of time scales between the elastic 



deformation due to flexure and that which results from the fluid's forcing (Tang 



and Paidoussis, 2008). The third parameter, k, is essentially the ratio of time scales 



between the elastic deformation due to fluid forcing and the elastic deformation 
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due to tension. These three parameters can be combined to produce: 



V 



U 



1/2 



KV 



Til 



(7) 



The utihty of which is that f //i effectively removes the dependence on sheet length 



from the scaled fluid velocity. This form was employed by Tang and Pai'doussis 



(2008) and it is also the appropriate velocity scaling for an infinitely long flexible 



plate under fluid loading (cf. Crighton and Oswell, 1991). The advantage of [3 



is that we have effectively removed the dependence on fluid velocity from the 
dimensionless tension so that we now have a direct measure of the deformation 
due to tension with that from flexure. This form appears in the study conducted 
by Yamaguchi et al.| ( |2003 ) for simply supported sheets. We should also point out 
that although /5 in ([t]) is derived by assuming that T acts in the axial direction, 
we can also employ this deflnition for (3 when T is applied at an angle 9 to the 
trailing edge of a three-dimensional sheet (cf. §[3]). In which case, T is understood 
to be the magnitude of the in-plane tension. 

We now construct a dispersion relation connecting the normalised circular fre- 
quency UJ and wave number k by placing ([2]) in (|5]), whence 



&{k,u)=uj^-^{k^ + (3e) + 2-l^^ 
f ^ \k\ 



[UJ 



kf = 0, 



(8) 



where j3 



KV 



A dispersion relation of similar meaning to (p|) is presented in 



Shelley et al. (2005) with axial tension omitted but intimated in the form of a 



Blasius friction law. By replacing f3 with f3 = 1.3R ^^'^L'^/B, where R denotes the 



Reynolds number, this later form appears in Connell and Yue (2007) for constant 



R. Here, however, we leave the magnitude of P unspecified. 
In solving the dispersion relation for uj we arrive at 



UJ 



± 



k 



\k\+2jj 



2/i± 



k'^\k\+ p \k\-2fiv'^ + 2fik'^ + 2fip 



1/2 



(9) 



for any real-valued k = nvr with n E N. It is interesting to note that ^ is not 
too dissimilar to the dispersion relation for a Kelvin-Helmholtz type instability of 
a two layer inviscid fluid, whereby the restorative influences of sheet flexure and 
axial tension are analogous to buoyancy and surface tension respectively. 

To determine instability, we search for frequencies which exhibit spatially grow- 
ing modes whereby ^{uj) > is required. Consequently, in ^ we set the terms 
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enclosed in square braces to zero and find that 



> (1 + |A;|/2/i) (e + p) 



(10) 



is required to realise instability of the flutter-type through convected flexural waves. 

Equations ^ and ( 10 ) with k = n are plotted in Figure [l| where we have 
chosen to scale the ordinate of each by fi which conveniently collapses the data. 
The role of f3 is to postpone instability; so that for f3 = 100, for instance, Vc/fi has 
increased by a factor of approximately 3.5 over (3 = 0. The fact that /5 does not 
affect u/fi at the point of instability will come as no surprise, given the form of 



For practical purposes we can write ( 10 ) in dimensional form. 



\K\ 
2pM 



iM + 2p/\K\){BK^ + T), 



which serves to illustrate the effects of both structural and added mass, flexure, 
and axial tension on stability. 

Despite their simplicity, equations (|9|, (10) and (11) illustrate some impor- 
tant aspects of the system. Firstly, any increase in the length of the sheet will 
destabilise the system — this effect is of 0{L^) which is quite disastrous from a 
practical standpoint. Both flexural rigidity and tension act to stabilise the sheet 
and any increase will be coupled by an increase in both Uc and u meaning the fluid 
must supply more energy to excite a fluid-elastic response. On the other hand, 
an increase in either structural or fluid added mass, determined by 2p/\K\, will 



promote instability. Another important point, first mentioned by Zhu and Peskin 



(2002) and more recently by Connell and Yue (2007), is that a massless system is 
infinitely stable. This is clearly demonstrated by (11) where f/c oo as M 



0. 



Lastly, it is worth pointing out that with axial tension replaced by a constant 



flow-induced tension from a Blasius boundary layer drag law (cf. Connell and Yue 



2007), it contributes P ~ C^(l) and therefore provides minimal additional stability 



when compared to an externally applied tension. 



3 Experiments 

We now outline the experimental campaign undertaken to investigate the drag and 
stability of a rectangular flexible sheet in uniform flow. The experimental campaign 
was conducted in the low speed wind tunnel of the Aerodynamics Laboratory at 
the Norwegian University of Science and Technology. The tunnel comprises an 
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available test-section of 2.7m x 1.83m and allows fluid velocities up to 30ms~^ with 
a turbulence intensity of less than 5%. The sheet was mounted vertically along the 
centre-line of the test-section with a clamped leading edge and free trailing edge. 
A diagram of the test set-up is provided in Figure |2] which we shall now explain in 
the following subsections. 

3.1 Apparatus 

A transparent Polyethylene material of thickness h = 0.15mm was chosen for the 
sheet as it provided a reasonably low mass per unit area m = 0.141kgm~^ and 
sufficient bending rigidity B = 59.65 x lO^^Nm^/m to allow an observation of 
flexural waves. Using a width of / = 75cm, the length of the sheet was varied 
to produce five aspect ratios from within the range 0.43 < ^ < 1. These five 
geometries are hereafter denoted SI, S2, S3, S4 and S5. The characteristics of 
which are summarised in Table [T] With these dimensions, Reynolds numbers were 
in the order of R = 10^ — 10^ — a range which encompasses the transitional regime 
for the boundary layer on a fiat plate. In addition to the aspect ratio, the other 
parameters varied were the magnitude of the in-plane tension T and the angle 6 
at which it was applied to the trailing edge. 

At the corners of the trailing edge, tension was applied. This was accomplished 
with two lengths of Kevlar string attached to two sets of known and equal weights. 
These weights were parametrically varied to produce a combined pre-tension in 
the order of 1 — 15N. To achieve a constant mean tension throughout a test run, 
each string travelled through a pulley system that directed it underneath the wind 
tunnel where it was attached to the weights. The position of the pulleys was varied 
to produce three different angles at which the string left the trailing edge: 6 = 0, 
22.5°, and 45°; whereby 6 = implies a perpendicular tension to the trailing edge. 
The tension in each string was measured directly behind the trailing edge of the 
sheet with two ring type strain gauges surrounded by swivels to avoid excessive 
rotations. The strain gauges were calibrated with an expected error in the order of 
iO.OlN. The applied tension produced a small oscillatory out-of-plane component 
when the sheet experienced flutter. At worst, its amplitude was around 7% of the 
total in-plane tension at f/ = 12ms~^ for S5. The magnitude of the out-of-plane 
component decreased for increasing such that for SI and S2 its amplitude was 
less than O.OIT at the maximum flow velocity considered here. 

A three-component force balance, situated underneath the centre of the test- 
section, was employed to measure the in-line force. This was accomplished by 
utilising a vertical pole, cantilevered from the force balance, that extended into 
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the centre of the test-section. Along the leading edge, nylon tape was employed, 
affixing the sheet to the pole and simulating a clamped boundary condition. To 
minimise its drag contribution and affect on the ffuid-elastic deformation of the 
sheet, the vertical pole, of height 1.43m, comprised an elliptical cross-section (a = 
4cm and b = 2cm). The bottom of each sheet was located 54cm above the base 
of the wind tunnel. This was deemed sufficient to avoid the developing turbulent 
boundary layer along the wind tunnel ffoor and walls — of thickness S ~ 15cm at 
the leading edge location, and, assuming a logarithmic velocity profile, climbing to 
approximately 17cm at the trailing edge position. 

The incident fiow velocity was obtained via a Pitot-static tube fixed 28cm from 
the roof of the wind tunnel. At an air temperature of 22.1°C, the properties of the 
fluid used throughout the experiments are p = 1.205kgm~^ and u = 0.150cm^s~^. 
Measurements from the Pitot-static tube, force balance and strain gauge channels 
were acquired through a HBM instrumentation amplifier attached to a notebook 
computer running customised data acquisition software. All data was digitised at 
200Hz and pre-processed through a low-pass Butterworth filter at 40Hz to avoid 
ambient noise. In addition, each run was recorded at 30 frames per second with a 
digital camera situated 44cm behind the trailing edge position on the wind tunnel 
fioor. 

Initial testing revealed the drag coefficient of the mounting pole to be Cd = 
0.59. Hence, the independent drag contribution of the pole ranged from 0.003N at 
0.5ms~^ quadratically climbing to a maximum of 1.59N at 12ms~^ — the largest 
fiuid velocity examined here. We expect that the addition of the attached sheet to 
the pole has the effect of a splitter plate and will therefore reduce the drag contri- 
bution of the pole further. Using a circular cylinder and splitter plate configuration 



as a guide (Hoerner, 1965, § 3-6), we can infer a conservative lower bound of O.QCd 
for the portion of the pole comprising the attached sheet. We can therefore expect 
a drag reduction by the presence of the sheet by, at most, 21% of the total pole 
drag. To place this into perspective, this amounts to, at worst, approximately 3% 
of the total measured drag force of the sheet at U = 12ms~^. Moreover, we suspect 
that given the elliptical cross-section of the pole, combined with the fiexibility of 
the sheet allowing some momentum exchange in the wake at the onset of fiutter, 
the drag reduction experienced by the pole will be somewhat smaller than our 21% 
upper bound. This effect is not corrected for in the data analysis implying that 
the calculated drag of the sheet will be slightly smaller than expected. 

A typical test run commenced from rest and the ffow velocity of the wind tunnel 
gradually increased over increments of approximately 0.2-0.5ms~^. This continued 
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beyond the onset of the first observed convected wave mode, until wild oscillations 
of the trailing edge ensued. At that point, the flow velocity was either terminated 
or gradually decreased over increments to cover the fluid-elastic instability region. 
Each run lasted for approximately 12-20 minutes, with time windows for each 
chosen fluid velocity, on average, 42.1s in length with a standard deviation of 
10.6s. This process was repeated over the parameter space of the test matrix which 
included five sheet geometries, three angles, and four to seven different magnitudes 
of trailing edge tension. 



3.2 Data analysis and processing 

Each data set was separated into time windows of length determined from the 
velocity profile incident on the sheet. After windowing, a general Fourier transform 
was performed to determine the average flow velocity, in-line force, and trailing 
edge tension for each time window. Moreover, the available frequency information 
allowed us to determine the dominant harmonics of the trailing edge tension and 
in-line force. The fluid-elastic instability was clearly observable from the trailing 
edge tension where the magnitude of the first harmonic increased by a factor of 
approximately 10 over one velocity increment once stability was lost. At this point, 
the Strouhal number, based on the length of the sheet (St = fL/U), was in the 
order of 0.4 — 0.6 — around half the value reported by Shelley et al.| (2005) for a 
cantilevered sheet in water with (3 = 0. As the flow velocity was increased beyond 
critical, St increased in a linear fashion to a peak of around 2.2 — 3.0 with the lower 
limit occupied by data sets comprising the largest in-plane tensions explored here. 

The drag experienced by the sheet was determined from the in-line force mea- 
sured at the force balance with the drag from the pole and total horizontal tension 
at the trailing edge subtracted. The drag force, denoted F^, was then normalised 
according to Co = Fc/O-bpAU'^ where A is the surface area of the sheet and U 
the flow velocity. The flow velocity, determined from the pressure measured by a 
Pitot-static tube, over each time window comprised a standard deviation of less 
than 0.05ms~^ from the mean for U > 2.4 ms~^. Below this range, however, the 
flow velocity showed some variability with the standard deviation increasing to 
0.19ms~^ in some cases. 



4 Results and discussion 

A selected run from the experiments is illustrated in Figure |3] which shows the 
incremental increase in fluid velocity and subsequent oscillatory behaviour of the 
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forces measured in-line and at the trailing edge once stability is lost to flutter. The 
behaviour illustrated in the plot once stability is lost, matches our observation of 
surface waves convected from the leading to the trailing edge of the sheet. The 
time series results presented in Figure [3] are quite typical of all cases considered 
throughout the experimental campaign. Frequency analysis of the time-windowed 
data indicates that the first modal frequency of the dominant convected wave event 
is in the order of 1-lOHz at the point of instability. The effect of which is transferred 
to both the force balance and strain gauges at the trailing edge (cf. Figure [s]). 
Rather than growing in time, these oscillations, which are steady and well defined, 
grow in amplitude and frequency as the fluid velocity is increased further beyond 
the first observed instability. When this occurs, a linear dispersion relation can 
no longer relate the frequency and wavenumber of the fiexural waves as elastic 
nonlinearities, such as curvature, become important. These nonlinearities will act 
to limit the amplitude of the instability. One could argue that this behaviour is in 
fact another stable state occupied before the sheet undergoes chaotic motion for 
flow velocities much larger than critical. However, we are not concerned with this 
chaotic motion at large fluid velocities here. 



4.1 Fluid-elastic instability 

We first examine loss of stability due to flutter. The critical velocity, denoted by 
f c, was determined by examining the first-harmonic of the trailing edge tension. A 
sharp jump in magnitude and phase shift of the first harmonic indicates stability 
loss and the appearance of convected waves. Figures |4],|5] and |6] illustrate the results 
obtained from this analysis for S2, S3 and S4 respectively. In each figure, the scaled 
velocity fc//U is plotted against the trailing edge tension for 6 = 0, 22.5°, and 
45°. We represent the tension by the parameter combination jS/fi"^ = T{M'^/ p^B), 
where T is the tension applied by the weight and the multiplicative factor M^/ p^B 
is constant throughout; the utility of which removes the dependence of sheet length 
and fluid velocity from the scaled tension and allows us to examine the effect of 
pure tension across each sheet considered. For comparative purposes, we include 



the stability boundary determined from the linear dispersion relation (10) for the 
first oscillatory mode k = tt. 

It is important to note that the exact instability boundary is not definitive in 
practice. There is a small but finite range of fluid velocities at which the sheet 
can undergo periods of both stable and unstable oscillatory motion. The trans- 
feral between these two states does not appear to be related to any of the input 
variables examined here. This bi-stability could be memory-dependent, such as 
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that demonstrated numerically by Zhu and Peskin (2002) and experimentally by 
both Watanabe et al. (2002) and Souilliez et al. (2006), or perhaps it may originate 
from perturbations of the flow velocity and any residual turbulence available. The 
dynamics of the system within this region are unclear and, in Figures |4| |5] and [6} 
we account for this uncertainty by incorporating a shaded region where upper and 
lower bounds of Vc are indicated by shape preserving polynomial interpolants. 

The experimental results for Vc indicate that instability can be postponed by 
increasing the in-plane tension. In most cases we observe a monotonic increase in 
fc//i with increasing Apart from a persistent vertical offset, this trend is 

partially captured by the qualitative expression for Vc (10). We suspect that this 
offset emerges from our limited description of the fluid added mass which appears 
in the dispersion relation. The qualitative theory assumes that the fluid pressure 
is composed of one single spatial mode that obeys only a periodicity condition 
at the sheet end-points. In reality, however, the the first modal frequency of the 
instability will comprise a number of spatial mode shapes. 

In contrast to the above, the results pertaining to S2 with tension applied at 
= (cf. Figure |4^) follow an entirely different trend. In particular, Vc measured 
at values of = and 400 are exceptionally large in comparison to their 

companion data. The dynamics of the instability appear to be different and we 
present Figure [s] to illustrate snap-shots of S2 with = 400 both before and 

after the appearance of convected waves. Although the left frame of Figure [8] shows 
S2 to be stable to convected waves, we observe a transverse oscillation mode. The 
right frame, at a higher fluid velocity, clearly shows this transverse standing wave 
envelope which is now superimposed on a convected wave instability. Its wave 
length is approximately 1/2. These well defined transverse modes were observed 
only for the shortest sheets, £^ = 0.75 and 1, where static deflection of the sheet due 
to its weight is minimal. Not surprisingly, with in-plane tension applied laterally 
{6 = 22.5° and 45°), these transverse modes are suppressed (cf. Figures |4]d and 
Souilliez et al. ( |2006 ), in the context of flags with /? = 0, suggests that 



such transverse modes and three-dimensionality can delay the onset of flutter by 
artificially increasing the rigidity of the system. This may well be the case because 
out-of-plane bending will increase the stiffness of the sheet to in-plane bending or, 
in our case, instability. 

With the magnitude of the in-plane tension held constant, an increase in 6 
postpones instability. To emphasise this effect, we present Figure [7] which shows a 
surface plot of Vc/ fi for S3 using a cubic polynomial interpolant over the measured 
data presented in Figure|5| The distinct trend of increasing stability with increasing 
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6, as demonstrated in the surface plot, is clearly reflected in the measured data 
for S2 and S3 (cf. Figures |4] and [s]) — note the values for fc at = and 

400 in Figure |4^ defy this trend, but, as previously mentioned, the dynamics of 
the instability are different here because a transverse mode shape was excited. 
Moreover, if one excludes the values of Vc corresponding to jS/fi^ = 2200 in Figure 

then S4 follows this trend also. For increasing 6, the in-plane tension comprises 
an increasing lateral component. Assuming that this suppresses three-dimensional 
geometric irregularities, or perturbations, then one can expect the sheet to adopt 
a more streamlined orientation and, consequently, an elevated stability. However, 
this trend is unlikely to continue for 6 > 45° where increasingly less axial tension 
would be available to the sheet. 

The effect of ^ on the sheet stability is difficult to determine from the present 
data set. In most cases, Vc is of the same order of magnitude for each sheet and 
no conclusions can be drawn. 



4.2 Fluid dynamic drag of a flag 



We now consider the fluid dynamic drag experienced by flexible sheet with (3 = 
— in other words, a flag. Measured results are plotted in Figure |9] for Cd = 
-Fo/O.SpAf/^ versus the Reynolds number, where the error bars of each data point 
correspond to the standard deviation of Cd for each time window. The results 
are presented for sheets with aspect ratios =2/ = 1, 0.75, 0.6, 0.5, 0.43. In addition 



to the present measurements, we compare with published data (Carruthers and 



Filippone, 2005), denoted 'CF' hereafter, for a cotton flag comprising aspect ratios 



= 0.1, 0.05, 0.03. 

For comparative purposes. Figure [9] includes the expected fluid drag of flat plate 



experiencing a turbulent boundary layer in plane parallel flow (cf. Schlichting and 



Gersten 2000, pages 580-586): 



CD 



InR 



G(A;D) 



with 



A = lnR, D = 21nK + K(C+ -3.0), 



(12) 



(13) 



where the function G(A; D) is determined from the numerical solution of (Schlicht- 



ing and Gersten [2000 , eq. 17.60) 

A 



, A 



D = A. 



(14) 
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This friction law is valid for smooth and rough surfaces by judicious choice of C+ 
for R > 10^. Given that the sheet material is polyethylene, it is reasonable to 
assume a smooth surface and therefore = 5.0. The constant k, denotes the 
Karman constant and takes the numerical value 0.41. 

As the aspect ratio £/ decreases, we observe a decrease in Cd (see Figure [9]). 
This trend is also apparent in the data of CF who studied ^ < 0.1 — a range 
more closely associated with ribbons and streamers. 

Despite the fact that our Co measurements are of the same order of magnitude 
and follow similar trends as those of CF, a slight vertical offset in the influence of 
^ on Cd is evident — cf. Figure [9] and the data sets corresponding to ^ = 0.43 
and CF's =2/ = 0.1 for example. Although the surface area of sheets considered 
by CF are (9(10) smaller than those examined here, the mass ratios fi are quite 
similar for both studies (/i = 5.9 — 10.2 for CF and fi = 6.4 — 15 here). However, 
given that CF employs cotton flags, the reduced velocity range is of (9(10^) higher 
than the present study. This suggests a much greater influence on the response 
of the sheet from elastic deformations. As such, this may account for the slight 
vertical offset between the two data sets. 

At small Reynolds numbers our measurements contain a degree of scatter which 
is generally reflected in the magnitude of the error bars. This uncertainty can be 
attributed to two factors: the small magnitude of the drag measured at small R; 
and the sensitivity of Cd to the sheet's orientation to the flow. For instance, as 
the flow velocity is increased the sheet assumes a number of stable non-oscillatory 
positions before aligning itself with the flow in a fully extended position. These 
changes in orientation affect the frontal and plane form area exposed to the flow 
and therefore, by consequence, the fluid drag. In other words, at low fluid velocities 
it is reasonable to expect that the memory or initial conditions of the system can 
influence the drag at subsequent fluid velocities before the sheet is fully extended. 
Once fully extended, we observe flutter. 

An example, however, where a fully extended sheet remains stable to convected 
waves is demonstrated for S2 where ^ = 0.75. Here, Cd exhibits two distinct 
regimes centred around a Reynolds number of R ~ 4 x 10^. This boundary corre- 
sponds to a sharp rise in Cd from 0.05 to 0.18 which approximates to a scaled flow 
velocity of f/yU ~ 30. At this flow velocity the sheet loses stability to flutter (cf. 
Figure |4^). A similar phenomenon is demonstrated by S4 with £?/ = 0.5, however, 
the rise in Cd is far less pronounced. For the remaining sheets, once fully extended 
in the flow, stability was lost immediately to flutter. 

The measured results indicate that the drag experienced by the flag is generally 
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of (9(10) greater than that of flat plate of equal dimensions. Given the sheet's 
behaviour at small R before reaching a fully extended orientation in the flow, 
together with its flutter response at large R, this result not surprising. Within our 
range of Reynolds numbers, (9(10^ — 10^), we can expect a turbulent boundary layer 
over the surface of the sheet. We notice that, like a flat plate, the drag coefficient 
decreases with increasing Reynolds number. However, the measurements indicate 
that this feature is perhaps more pronounced than that experienced by a flat plate. 
This would suggest that the boundary layer thickness is larger for a flexible sheet. 



4.3 Fluid dynamic drag with in-plane tension 

We now focus our attention on sheets S2, S3, and S4 and examine the fluid drag ex- 
perienced by each in the presence of in-plane tension. We represent the magnitude 
of this tension by /? = T{L'^/B) which is applied at an angle 9 to the trailing edge. 
For S2 and S3, the drag coefficient is presented for = and 9 = 45° in Figures 



10 and 11 respectively. For S4, however, we consider an additional intermediate 



angle of ^ = 22.5° and these data sets are shown in Figure 12 For comparative 



purposes, the drag experienced by a flat plate (12) is plotted on each set of axes. 

The measurements demonstrate a significant drag reduction through a moder- 
ate application of in-plane tension at the trailing edge. For example, S2 with (3 > 
experiences a drag reduction by a factor of approximately Co/<i when compared 
to that of P = (cf. Figure 10). Moreover, such a dramatic reduction depends 
not only on the magnitude of the tension, but also on 9. To complicate matters 
further, stability also plays a crucial role. Distinctly different drag characteristics 
are observed before and after stability is lost to convected waves. 

Before instability v < Vc, the magnitude of Cd appears to be independent of 
(3 > 0. Provided sufficient in-plane tension is applied, the fluid drag of the sheet is 



reduced regardless of [3 (see Figure 11 d in particular). However, as 9 increases the 
drag experienced by the sheet is further reduced. This behaviour is presumably 
caused by an improved and more streamlined orientation of the sheet — this relates 
to the discussion involving an increased fluid-elastic stability for increasing 9 (cf. 



4.1 ). This effect is best demonstrated in Figure 12 for S4. At = 0, the measured 



data is perhaps twice the Cjj expected for a flat plate. While for 9 = 22.5° and 45°, 
the drag is reduced to such an extent that it is given to reasonable approximation 
by the flat plate drag law. 

In contrast, for f > fc, we observe that the in-plane tension plays a crucial role 
in the drag characteristics of the sheet. Ultimately, the magnitude of the in-plane 
tension dictates the critical fluid velocity, and the subsequent rise in the fluid drag 
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once stability is lost. (cf. Figure 10 d). This explains the cascading behaviour of 
Coi for increasing (3 during increasing Reynolds numbers, demonstrated in Figures 
[T0| 3 and [TT)d. Furthermore, these figures also suggest that the drag behaves in a 
linear fashion once stability is lost. In particular, we observe a six fold escalation in 
Cd past the critical velocity for a moderate Reynolds number increase of 10^. The 
convected wave response of the sheet is responsible for this dramatic drag increase 
once the critical velocity is breached. Although not indicated over the range of 
our data, one would expect that this linear increase in Co to be arrested once 
the sheet performs chaotic motions and the coherent structures in the wake of the 
sheet breakdown (see Connell and Yue, 2007, Fig. 23 and the related discussion) 
— this was not investigated here. 



5 Conclusions 

In the present work we have experimentally studied the drag and stability of a 
cantilevered flexible sheet in uniform flow characterised by R = 10^ — 10®. For the 
canonical case of a flag, we observe a decrease in the drag coefficient for decreasing 
aspect ratios within the range 0.43 < l/L < 1 — this agrees with previous studies. 
In the context of fluid drag and fluid-elastic stability, the notable feature of this 
present work has been an investigation into the effect of in-plane tension applied 
to the sheet. Our results indicate that the onset of a flutter- type instability can be 
delayed through the application of an in-plane tension. In particular, a monotonic 
increase in the critical velocity is observed for increasing tension. In the presence 
of in-plane tension, a dramatic reduction in the drag coefficient of the sheet is 
observed. The unique drag characteristics of the sheet are a corollary of its fluid- 
elastic stability. Interestingly, two regimes in the drag coefficient are seen, which 
appear to be separated by the critical velocity. Before stability is lost, the drag 
is well behaved and, provided transverse modes are suppressed, given to good 
approximation by the turbulent boundary layer drag law for a flat plate. However, 
once the sheet loses stability, the drag coefficient increases rapidly in a linear 
fashion with respect to the onset flow. The point at which this rapid increase 
occurs depends on the critical velocity and therefore the magnitude of the in-plane 
tension applied to the sheet. 
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Figure 1: Left axis: the critical velocity Uc///. = U{M^ / p^By/'^ plotted against the 
mass ratio n — pL / M ior the first eigenmode k — tt oi a flexible sheet under axial 



tension P 



KV 



in uniform flow. Right axis: the corresponding circular frequency 



oj/ jji = Vt{M/ pU) at the point of instability. 
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Figure 2: A schematic of the test set-up. The wind tunnel floor is represented by 
the grey horizontal line. 
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Figure 3: A typical time series for S2 
22.5°to the trailing edge, (a) illustrates 
the forces measured: in-line Fx; at the 
strain gauge T2. 
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Figure 4: S2: £?/ = 0.75. The measured critical flow velocity fc//i = t/(M^/p^i?)^/^ 
plotted against trailing edge tension = Ti^M"^ / p^B) applied at angles of: (a) 
^ = 0; (b) ^ = 22.5°; and (c) 9 = 45°. The dashed line corresponds to Eq. (|To|. 
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Figure 5: S3: = 0.60. The measured critical flow velocity fc//i = ?7(M^/p^i?)^/^ 
plotted against trailing edge tension = T{M'^/ p^B) applied at angles of: (a) 
^ = 0; (b) ^ = 22.5°; and (c) 6 = 45°. The dashed line corresponds to Eq. ([To]). 
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Figure 6: S4: = 0.50. The measured critical flow velocity fc//U = f/(M^/p^i?)^/^ 
plotted against trailing edge tension (3/ p'^ = T{M'^/ p^B) applied at angles of: (a) 



^ = 0; (b) ^ = 22.5°; and (c) 9 = 45°. The dashed line corresponds to Eq. (10). 
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Figure 7: A surface plot of the measured critical flow velocity v^/ IJ. — 
U{M^/p'^By/^ versus the in-plane tension P/fj,^ = T{M^/p'^B) and its direction 9. 
Cubic polynomial interpolates have been used along the P/p"^ axis to re-grid the 
measured data for mean values of Vc- This plot pertains to S3 where S2f = 0.6. 



24 



Figure 8: S2 with trailing edge tension (3 = 400 applied at 6 = 0. A transverse 
mode is observed across the surface of the sheet. Left frame: v < Vc] and the Right 
frame: v > Vc- 
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Figure 9: The fluid drag in uniform flow experienced by a horizontal flag. The 
present measurements for a polyethylene sheet are plotted for five aspect ratios 
^ = l/L with error bars corresponding to the standard deviation of Cu at each 



measured flow velocity. Published data (Carruthers and Filippone, 2005), denoted 
'CF', for a cotton flag provide three additional values of . The solid line rep- 
resents the expected fluid drag of a flat plate experiencing a turbulent boundary 
layer (cf. Schlichting and Gersten, 2000, eqn. 18.99). 
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Figure 10: The measured drag of S2 under constant tension (3 = kv'^ (x 10^) in 
uniform flow. The subplots correspond to the trailing edge tension applied at 
angles: (a) 6 = 0; and (b) 6 = 45°. The error bars for each data point correspond 
to the standard deviation of Cd at each measured flow velocity. The solid line 
represents the drag on a flat plate, of j^/ = 0.75, experiencing a turbulent boundary 
layer. 
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Figure 11: The measured drag of S3 under constant tension j3 = kv'^ (x 10^) in 
uniform flow. The subplots correspond to the trailing edge tension applied at 
angles: (a) 6 = 0; and (b) 6 = 45°. The error bars for each data point correspond 
to the standard deviation of Cd at each measured flow velocity. The solid line 
represents the drag on a flat plate, of ^ = 0.6, experiencing a turbulent boundary 
layer. 
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Figure 12: The measured drag of S4 under constant tension (3 = kv"^ (x 10^) in 
uniform flow. The subplots correspond to the trailing edge tension applied at 
angles: (a) 6 = 0; (h) 6 = 22.5°; and (c) 6 = 45°. The error bars for each data 
point correspond to the standard deviation of Cd at each measured flow velocity. 
The solid line represents the drag on a flat plate, of ^ = 0.5, experiencing a 
turbulent boundary layer. 
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Table 1: Characteristics of the sheets employed in the test matrix where denotes 
the aspect ratio, A the projected surface = pL/M the mass ratio, v — 

U{ML?' /By/"^ the reduced or normalised incident flow velocity, kv"^ = TL^/B the 
ratio of tension to flexure, R the Reynolds number, and 6 represents the angle of 
application of tension at the trailing edge. 



Sheet 




A (m^) 


1^ 


V 


KV^ ( 


xlO^) 


R (x 


105) 


e 


SI 


1.0 


0.563 


6.410 


23.0 - 


434 


19.4 


- 75.0 


0.315 - 


- 5.95 





S2 


0.75 


0.750 


8.546 


37.4 - 


549 


34.5 


- 133 


0.513 - 


- 7.53 


0, 22.5°, 45° 


S3 


0.60 


0.938 


10.68 


43.2 - 


599 


54.0 


- 208 


0.592 - 


- 8.21 


0, 22.5°, 45° 


S4 


0.50 


1.125 


12.82 


56.2 - 


744 


77.7 


- 300 


0.770 - 


- 10.2 


0, 22.5°, 45° 


S5 


0.43 


1.312 


14.96 


63.8 - 


977 


106 


-408 


0.875 - 


- 13.4 
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